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Abstract 

We show that some models with non-local (and non-localizable) interactions 
have a property, called quasi-locality, which allows for the definition of a transfer 
matrix. We give the Yang-Baxter equation as a sufficient condition for the existence 
of a family of commuting transfer matrices and solve them for a loop model with in- 
tersections. This solvable model is then analyzed in some detail and its applications 
to a Lorentz gas are briefly discussed. 

In the field of solvable models in statistical mechanics it is customary to restrict oneself to 
models with local interactions. This restriction makes it hard to treat models of extended 
objects, such as linear or branched polymers, directly. The usual approach is to consider 
the graphical expansion of a bona-fide model with local interactions, such as the Potts 
model or an 0(n) model, and derive properties of the graphical model from the local 
model. This is not always possible. In this paper we present a model, the intersecting 
loop model, for which there is no local representation as far as we know. 

We consider a dense gas of loops on a square lattice. In contrast to the loop gas that 
one obtains in the graphical expansion of the Potts model or a RSOS model, our loops 
are allowed to intersect themselves and each other. 

On every site of a square lattice we randomly put a mirror along one of the diagonals, 
or leave it empty. The mirrors deflect photons that run along the edges of the lattice. If 
we impose toroidal boundary conditions every photon follows a closed orbit. Obviously 
every edge belongs to precisely one path. In this paper we are mainly interested in the 
statistics of the paths, i.e. given the probabilities that a site is empty or that it carries a 
mirror oriented along the NE-SW or the NW-SE diagonal we want to calculate quantities 
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such as the density of paths and their length distribution. The kinetic properties of this 
Lorentz gas, such as its anomalous diffusive behaviour, have been sudied in |l] . 

We take the probabilities P^NE-SW), Pj(NW-SE) and P { (empty) that the mirror at 
site i is oriented along the NE-SW diagonal, the NW-SE diagonal or absent independent 
of i. For convenience we will mostly work with the Boltzmann weights (unnormalized 
probabilities), denoted respectively by a, b and c. Furthermore we give every loop a 
fugacity q. The oject of our study is the partition sum (or generating function) 

% — a #NE-SW & #NW-SE c #empty^#paths_ ^ 

mirror 
configurations 

When c = this model can be mapped onto a 6 vertex model and it is the graphical 
representation of the self-dual g 2 -states Potts model. However, if c 7^ no mapping to 
a model with local interactions is known, unless q is an integer > 1, hence a different 
approach is needed. 

Our approach is based on the observation that, although the Boltzmann weight is 
non-local, it is possible to obtain the Boltzmann weight for a configuration on a large 
patch of the lattice by glueing together smaller patches. This property may be called 
quasi-locality. 

Let A and B be neighbouring patches of the lattice. To every mirror configuration 
on A we assign a Boltzmann weight VT(A) (a factor a, h or c for each site and a factor 
q for each closed loop in A) and a connectivity, indicating which pairs of edges on the 
boundary dA of A lie on the same photon path. We do the same for B. Now we can 
calculate the Boltzmann weight and connectivity of the total patch AB. The connectivity 
of dAB is induced by the connectivities of dA and dB (glue the paths together on the 
common boundary dA n dB of A and B) , and the Boltzmann weight is 

W(AB) = W(A)W(B)q* {9AndB \ (2) 

where j^{dA R dB) is the number of paths closed by glueing A and B together. 

This property makes it possible to define a transfer 'matrix', acting in the space of 
connectivities of the boundaries of a cylinder. It is defined via 

Z M+l (a) =^2T aP Z M {f3). (3) 



Here a and (3 denote connectivities of the boundaries of a cylinder (the circumference iV is 
implicit in the equations) and Zm{®) is defined as the sum over all mirror configurations 
that lead to the connectivity a of the Boltzmann weight for mirrors and closed loops on 
a cylinder of height M. This definition of a transfer matrix is similar to the one for the 
Q-states Potts model in [0. 
Since 

Z (a) = 5 aid , (4) 

where id is the trivial connectivity, in which every edge in the top row is connected to the 
edge in the bottom row straight below it, 

Z M (a) = (5) 
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The toroidal partition function is obtained by closing the cylinder. This yields a factor q 
for every newly formed loop; for each connectivity a we denote the number of loops by 
#(«). Finally we sum over a: 

Z M = Y.^ (a) T^. (6) 

a 

We now rephrase these observations in a slightly more formal language. The con- 
nectivities of the upper and lower boundaries can be taken as the basis of an algebra. 
Addition is defined formally and multiplication is defined by glueing the connectivities 
together, again giving a factor q for each loop that is formed. Since there are (2iV — 1)!! 
connectivities, the algebra of connectivities B^{q) (a.k.a. the Brauer algebra, since it was 
discovered by Brauer |3j as the algebra of invariants of SO(q) for integer q) is (2iV — 1)!! 
dimensional. 

The connectivities are most naturally represented by pictures. Each picture consists 
of two rows of iV dots each above each other, representing the edges on the upper and 
lower boundary, and lines connecting pairs of dots, indicating which edges lie on the same 
path. To calculate the product a/3, identify the bottom row of a with the top row of /3, 
remove the dots on this middle row, straighten out the lines and replace every closed loop 
by a factor q. 

Instead of considering a connectivity and its partition function separately, we can view 
them as an element of the Brauer algebra: 

Z M {a) >-> Z M (a) ■ a. (7) 

If we now define 

T M = ^Z M (a) ■ a (8) 

a 

the transfer matrix, defined in eq. (0), becomes an element of the Brauer algebra, 

Z M +i = TZ M - (9) 
The relation with the 'components' T a p given before is 

Tp = J2T a p-a. (10) 

a 

We denote a configuration of a row of mirrors by fa = (mi, m.2, • ■ ■ , tun), where rrii = 1, 
— 1 or according to whether site i has a NE-SW mirror, a NW-SE mirror or is empty. 
To every fh corresponds a connectivity a(rh), so the transfer matrix can be written as 

T = y^to(m) • a(rh), (11) 

m 

where 

W (rn) = ? <5^ > 0) a #(m i =l) 6 #(m i =-l) c #(m i =0)_ (12) 

The transfer matrix is the linear combination of connectivities 



T = -(*)-(*) (*M«> 



(13) 



where the leftmost horizontal edge should be connected to the rightmost horizontal edge 
and where 



— (^y~ = a — ^/ — + b — ^ — + c — | — 



(14) 



Since Z n = id, 




Z M = T M . (15) 
The partition function Zm can now be written as 

Z m = TtZ m , (16) 

where the trace is defined as a linear function B^{q) — > (D with 

Tra = g #(a) . (17) 

The trace satisfies 

Tr(a/3) = Tr(/3a). (18) 

A simple example may clarify the definitions. Consider the algebra B±(q). It is 105- 
dimensional and two connectivities are 

W 1/ 

a= and (3= (19) 

Their products are 

o>(3 = q and (3a = q ~ _ ~ (20) 
and the traces are 

Tra = g 2 Tr p = q Tr(a/3) = Tr(/?a) = g 3 . (21) 

The usual lattice models with local interactions are called solvable when their Boltz- 
mann weights satisfy the Yang-Baxter equation. This equation implies that the transfer 
matrices constructed with the aid of these Boltzmann weights commute for arbitrary sys- 
tem sizes. We can do something similar for the loop model: when the mirror probabilities 
a, b and c satisfy the Yang-Baxter equation for loop models, the corresponding trans- 
fer matrices T commute for arbitrary system sizes. The Yang-Baxter equation for loop 
models states that the l.h.s. and the r.h.s. of eq. (p2|) are the same linear combination of 
connectivities of the six outer edges. 




Here 



X 



X*'X 



(23) 



Explicitly, the Yang-Baxter equation reduces to the following set of four equations 



cb'a" = bb'c" + bc'a" 

aca = aa c + ca a 

ab'c" = ac'b" + cb'b" 

ab'a" = qba'b" + aa'b" + ca'b" + bb'b" + bc'b" + ba'a" + ba'c" 



(24) 



Regarding these equations as a linear system for the unknowns a", b", c" and requiring 
the determinant of the coefficients of the first three equations to vanish leads to 



ab(a' + b')c' = a'b'(a + b)c. 
If we want solutions with positive a, b and c we must have 

ab a'b' 

K. 



c(a + b) d{a' + b') 
The fourth equation from the set fl24]) then implies that 

re = 1 -q/2. 

We can thus parametrize the mirror weights as follows 

a(u) = p{u) ■ (1 — u) 

b{u) = p(u) ■ u 

c{u) = p(u) • ku(1 — u), 



(25) 



(26) 



(27) 



(28) 



with re given by eq. (]2"T|). The function p(u) is an arbitrary normalization. Furthermore we 
will assume that < re < 1 so that the model allows for intersections and has a physical 
regime < u < 1 where a, b and c are positive. In the sequel we will take p(u) = 1, unless 
otherwise stated. The spectral parameters satisfy the difference property u" —u' — u. 

The proof that a solution of the Yang-Baxter equations leads to commuting transfer 
matrices, i.e. 

T(u)T(t;) = T»T(u) Vw,y (29) 

is the same as usual. It makes use of the periodic boundary conditions, the existence of 
a 'horizontal' inverse of the 'R-matrix' for generic values of the spectral parameter: 





(30) 



and repeatedly uses the YBE: 



T'T 




= TT'. 



(31) 



There are four points where the R-matrix R{u") doesn't have a horizontal inverse, viz. 
u" = ±1 and u" = ±k _1 . At u" = 1 and u" = k~ x it is a multiple of a projector, i.e. at 
these points there is a non-singular constant A such that the R-matrix satifies 



so that 



R ■ R — XR, 
X^R- (\6 2 - X^R) = 0, 



(32) 
(33) 



where \di2 denotes the identity element in the 'horizontal' Brauer algebra B-2{q). This 
implies that for these two values of u" the matrix products T(u + u")T(u) split in the sum 
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of two terms, one with local building blocks 




the other with building blocks 




For u" = 1 (and hence u' = u + 1) the first term ([34] ) enormously simplifies. At this point 
A = q and the YBE becomes 




Since on the r.h.s. the two points on the left are always connected, the same must hold 
for the l.h.s. It now follows that 



It is easy to see that 
Furthermore, 
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and the second set of configurations (35) turns out to have an overall factor u. We thus 
obtain the inversion identity 

J(u + l)T(u) = <j)(u) N \d + u N TV>(u), (40) 

where the 'fused' transfer matrix T^ 2 '(u) is a polynomial of degree 3iV in u over the Brauer 
algebra Bj^(q). 

As for many other models, for large N and small u the first term on the r.h.s. of 
eq. (0) dominates the second, hence every eigenvalue A(iz) of T satisfies 

^{u + l)V>(u) = <f>(u), (41) 

where 

Mu) = lim (A(u)) 1/N . (42) 

N— >oo 

In particular the 'partition function per site' ipo{u) satisfies the inversion relation (f|l]). It 
is fixed by the symmetry u —>■ 1 — u (related to an involutive automorphism of the Brauer 
algebra) to be a product 

Mu) =F(u)F(l-u) (43) 

and the function F{u) satisfies 

F(l + u)F(u) = + + (44) 
The solution that is free of zeroes in the 'physical' strip < Rew < 1 is 

Y(\ + —) T(- + — + -) 

F( U ) = 2^ " { ( ;ii ^I 2 ;V ■ (45) 

Hence the result for the normalized partition function per site (with a + b + c = 1) is 

, 4„ r(i + f)r(|-f)r(i + ^ + f)r(i + ^-f) 

^ 0n i + «u(i-«)r(i + s)r(i-f) r(i + f) r(i + £-f)- 1 j 

When we take the logarithmic derivative of the normalized partition function with respect 
to k, keeping u fixed, we get the linear combination of densities of mirrors and loops: 

9 \ na ,i, da ^ j_ dh ^ j- dc ^ j_ a 5 (Pioops)_ (A7 , 

^-lOg^On— o 1 - "o r~ + « • ^ 4 'J 

OK ok a OK OK c ok q 

When q — 1 (k — |) the densities of the various types of mirrors are are simply equal to 
their normalized Boltzmann weights, since there is no coupling between different mirrors 
(the partition function for q = lisZ = (a + b + c )# sites )_ Then eq. ([|7]) reduces to 



(Ploops)g=l = ~ ^log^On 



(4* 

K=l/2 



At the isotropic point u = |, which corresponds to the mirror weights 



, 4 1 

6 =9- C =9' (49) 



the loop density is 

(pioo PS ) g =i = 2 (V 2 (f) -Ml)) -| 

= 27r-f, (50) 

with ip2 the digamma function. This number can be measured in a Monte Carlo simula- 
tion, by setting the mirrors according to the probability distribution eq. ([49]) and counting 
the number of paths. The numerical result agrees with the exact answer to five significant 
digits 1. 

The loop density also gives information about the length distribution of loops passing 
through a given edge: the average inverse length is given by 

(L- 1 ) = i(pioo P s) 9 =i. (51) 

A more extensive and detailed study of this loop model will be published elsewhere [|5| . 
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